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Abstract: 



A magnetic susceptibility which decreases with decreasing temperature is observed in all Cu02 
based superconductors with less than optimal doping. We propose that in La2-xSrxCuOi this is 
due to antiferromagnetic ordering which is prevented by the low spatial dimensionality while in 
YBa2Cu30Q,Q it is due to the interplay between antiferromagnetic fluctuations within a plane and 
singlet pairing of electrons between nearest neighbor planes. 
Pacs: 74.65+n,76.60-k,75.30.Kz,75.10-b 
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The low frequency spin dynamics of the CM02-based superconductors are anomalous [^. The 
anomaly which has received the most attention is the different temperature dependences of the 
copper (Cu) and oxygen (O) nuclear magnetic relaxation rates Tf^. The difference apparently 
occurs in all CUO2 superconductors and has been successfully interpreted |@] in terms of antifer- 
romagnetic spin fluctuations which increase in strength as the temperature is lowered, and which 
couple much more strongly to the Cu than the O nuclear moments. It is still controversial [|^, ^, ^, |^ 
whether the increase in strength is associated with a temperature dependent correlation length as 
expected near an antiferromagnetic ordering transition. A less well understood anomaly occurs 
in "underdoped" CUO2 materials (i.e. those with less than the optimum number of carriers) and 
involves a decrease in the spin susceptibility Xs and some relaxation rates 1/TiT with decreasing 
temperature, as expected in systems with a gap to low-lying spin excitations. Several authors have 
proposed that this spin gap behavior is evidence for the existence of exotic quantum liquid phases 
in underdoped high materials P]. Most of the theoretical work on exotic quantum phases has 
focussed on a possible doping induced spin-gap state of a single CUO2 plane. In this paper we 
argue that the data suggest a different physics: in La2-xSrj.Cu04 at x < 0.15 an antiferromagnetic 
instability which is prevented from developing into true long range order by the low dimensional- 
ity nevertheless leads to formation of a gap in the longitudinal spin fluctuation spectrum which 
leads to a decrease (but not to zero) of Xs- However, in YBa2Cu30Qs the physics is driven by 
an interplay between tendencies towards spin singlet pairing of electrons in adjacent CUO2 planes 
and antiferromagnetism in a given CUO2 plane. 

The ideas we present here are related to previous work on the high problem. D. C. Johnston 
argued early on that a Heisenberg model with a J decreasing with increasing x explained the 
La2-xSrxCuOi data 0]. Magnetic correlations are also important in the spin bag picture of 
Schrieffer, Wen and Zhang Barzykin and Gorkov have discussed theoretically the possiblity 
of strong short-ranged magnetic correlations in La2-xSr.j.CuOi and mentioned the connection of 
this with a Xs{T)- Our analysis of Y Ba2Cu-iOQ,Q is related to previous work of Altshuler and loffe 
as discussed below |10[ . 



We now turn to the data. Xs may be determined from the measured bulk susceptibility 
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Xb = Xs + Xc + Xvv + Xdia if the core (xc), van Vleck (xw) and Landau diamagnetic 
iXdia) contributions are known. In YBa2Cu307, Xdia has been found to be less than 1/6 of Xs 
[pT| . We neglect Xdia in our analysis. Xs for the Y Ba2Cu^0-j-s materials is known [|l], For 



La2-xSrxCu04, however, there are conflicting assertions in the literature |]I3|, |TJ, |1^, |l6l. We 
estimate that in La2Cu04^ {xc + Xw) = l-Ostates/{eV — Cu) for fields parallel to the Cu02 
plane and —l.?)States/{eV — Cu) for fields perpendicular to it by comparing the measured [0 
susceptibility (at temperatures of order 600-800K, well above the 3d Neel ordering temperature) 
to the known values for the 2d S=l/2 quantum Heisenberg model [|T^. The example of the 



Y Ba2Cu^0T-s series |rT| and the near equality [|13| of the difference between the ab-plane and 
c-axis susceptiblities of La2CuOi and Lai,ij2bSrQ_Q-j^CuOi shows that it is reasonable to assume 
that Xvv and Xc do not change much with doping; we have used this assumption to extract Xs from 
Xb- 

Fig. 1 shows Xs{T) for some high materials, and for the S=l/2 2d Heisenberg model with 



J = 0.13 eV. Fig. 2 shows the copper nuclear relaxation rate (1/TiT) data [T^, ^ It is clear that 
YBa2Cu30Q_Q is different from the other materials. Its spin susceptibility is smaller in magnitude 
and larger in relative temperature dependence than all of the others and unlike the other materials 
extrapolates to zero at T = 0. The relaxation rates in La2-xSrxCu04^ are larger and more 
strongly temperature and doping dependent than in YBa2Cu-^07^s- Further, in La2-xSrxCuOi 
the Cu 1/TiT increases as T decreases except in a small region below about T=50K (which we 
suspect is dominated by superconducting fluctuations), in contrast to Y Ba2CuzOQ,Q where the Cu 
1/TiT has a broad maximum at T=150K. In YBa2Cu30Q_Q the oxygen 1/TiT has nearly the same 
T-dependence as XsiT) [||, [III] . 



The distinctive behavior of La2-xSrxCuOi, namely a decrease, but not to zero of Xs and a 
monotonic increase in 1/TiT as T —>■ 0, are signatures of antiferromagnetism; e.g. Xs for the 
2D S=l/2 Heisenberg model is maximal at a temperature of order J and drops by about a 



factor of two between T ~ J and T = |T^. However the behavior of Lai,s6Sro,i4Cu04 



cannot be interpreted in this way: although Xs{T) is at least qualitatively consistent with that 
of the 2D S=l/2 Heisenberg model with J of ~ 300 - 400K, the 1/T behavior of the Heisenberg 



model at T > J is not observed at T > 400-?^ [|T3[. A "two- fluid" model [0 with only a small 
density of Heisenberg spins is not compatible with the relaxation rate data which imply rapid 
T-dependent relaxation of all Cu spins. However, fluctuating regions of local spin density wave 
(SDW) ordering can explain the data. To see this, note mean fleld theory predicts a transition at a 
temperature Tmf to a phase with nonzero staggered magnetization N. In two spatial dimensions, 



thermal fluctuations prevent long range order for T > for T < Tmf the appropriate 
picture is of slowly fluctuating domains, with non-zero, but random from domain to domain. 
We believe Xs in this situation may be reasonably approximated by rotationally averaging the 
mean fleld result (which depends on the angle between the fleld and N), and rounding out the 
singularity at Tmf- This leads to a susceptiblity which drops by a factor of 2/3 between Tmf 
and 0. This calculation neglects quantum effects which would reduce Xs further. These will 
be less important as Tmf is decreased, because a lower Tmf implies a larger bare correlation 
length at Tmf and hence a large "effective spin" which is ordering. These arguments imply that 
antiferromagnetic fluctuations at T < Tmf lead to a Xs(T). However, a mode-coupling analysis of 



antiferromagnetic fluctuations at T > Tmf yields a negligible T-dependence of Xs f^^- Thus we 
propose that in the Lai,g25Sro,Qi5Cu04 sample the Tmf is rather above room temperature while in 
the Lai sQSrQ i4^Cu04 sample it is somewhat below. This may be consistent with neutron scattering 
experiments on Lai g^Sro 14CUO4 |Q in which quite sharp peaks are observed emerging at low T. 
For still larger Sr concentrations, Tmf < T^ and Xs displays a very weak T dependence. 

The Cu NMR relaxation is dominated by the relaxational dynamics of the AF fluctuations, 
so 1/TiT ~ J2g^^^^ fi^O/oi^lO ~ where ^ is the correlation length, 77 and z are critical 

exponents and f and g are scaling functions for the staggered susceptibhty and spin fluctuation 
energy scale respectively. Note this form is not multiplied by Xs{T) in contrast to ref ||16|. The 
"MMP" form proposed for yi?a2CM307 corresponds to ?7 = z = 2, f = = + 
For T > Tmf, the MMP forms are appropriate; as T ^ the function must cross over to the 
2D AF scaling forms, where 77 = 2 and z = 1. The oxygen relaxation rate, however, is due to 
small q spin fluctuations, i.e. to fluctuations of the (nearly) uniform magnetization Mg In an 
ordered antiferromagnet relaxation is due to spin waves; At low T the number of thermally excited 
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spin waves is small and the projection of these onto Mq vanishes as g — > 0, so that the oxygen 
1/TiT ~ [^. Similarly in a 2d SDW below Tmf the formation of antiferromagnetic domains 
will lead to an oxygen 1/TiT which drops more rapidly than Xs{T)i in contrast to the prediction 
in [0. If our analysis of XsiT) is accepted then this behavior has already been observed [T^, 



Now consider the Y Ba2Cu^0'j^s system, where the basic structural unit is a pair of Cu02 
planes, separated from the next pair of Cu02 planes by the relatively inert CuO chains |^6[. We 



model the spin dynamics of this system using two coupled planes of antiferromagnetically correlated 
spins. This is a gross oversimplification of the physics of the real material, because it omits the 
charge degrees of freedom. We discuss this further below. The Hamiltonian is 

<i,j>,a i 

Here i and j label nearest neighbor sites in a two dimensional square lattice and a G {1,2}, 
labels the two different planes. This model has a T = phase transition between a large J2 singlet 
state and a small J2 antiferromagnetic state. To study the regime near the transition we use the 
Schwinger boson method . In the model of Eq. |l] we find the second order transition of interest 
here to be preempted by a first order transition. For our detailed calculations we used a slight 
variant of Eq. |I| in which the first order transition is suppressed. In the mean-field analysis, sums 
over the momentum q occur; these may be replaced by an integral over an energy times a density of 
states which, for the model of Eq. 1, is constant near the band edges and logarithmically divergent 
at band center. Replacing this density of states by a constant yields a model with a second order 
T=0 transition (of the 3D Heisenberg universality class) at = 4.48Ji. Insulating YBa2Cu30QX) 
presumably has Ji ^ J2 [^. We believe that as the doping increases, the effective Ji decreases 
and J2 increases. In YBa2Cu30Qs antiferromagnetic correlations between pairs of planes have 
been observed to persist up to room temperature . 

We have computed the temperature dependence of {xs{T)) and of the oxygen, yttrium and 
copper NMR relaxation rates for various J2 > J2 using Eq. 0, the constant density of states, and 
a simplified version of the standard NMR form factors in which the Cu transferred hyperfine 
coupling B was set to 0. Some results are shown in Fig. 3 for J2 — J2 = 0.3 ; the resemblance 



of the curves for Cu and O to the data for YBa2CusOe,Q is evident. Note in particular existence 
of two temperature scales; a higher one, of order J2, at which Xs and the oxygen 1/TiT begin to 
drop, and a lower scale of order (J2 — J*) < J2 at which the Cu 1/TiT begins to drop. The 
different T dependences of the Cu and O 1/TiT in our model are due to both a growing correlation 



length and different spin gaps in different regions of q-space, in contrast to the previous model 
in which the effects of the spin gap were modelled by multiplying the MMP x by Xs{T), and in 
which the difference was due only to a T dependent correlation length. The difference between the 
O and Y relaxation rates arises in our model because the Y nucleus is relaxed only by fluctuations 
symmetric under interchange of the two planes; these are the most strongly suppressed by the 
tendency to form singlets. The decrease of the yttrium 1/TiT relative to Xs comes from the phase 
space restriction involved in scattering from one small q low-energy state to another and may be 
a special feature of this spin-only model. The difference between the O and Y relaxation rates is 



not consistent with published Y relaxation rate data ||3^ but the uncertainty in the Y relaxation 
rate is large. 

We have also computed the static susceptibitlity at g = {n,n), XApiT), finding it to be only 
very weakly T-dependent for J2 — J* = 0.3, because of the competition between singlet formation 
(leading to a xaf decreasing with decreasing T) and antiferromagnetism (leading to the opposite). 
For all J2 > J2 we find xaf decreases with T for T less than the temperature at which 1/TiT for 
Cu has its maximum. This follows from a Kramers-Kronig argument: a decrease in 1/TiT implies 
a shift of spectral weight in x"(g,ti;) from lower to higher uj; x'(g,u; = 0) = / dujx"{q,uj)/uj must 
then decrease. It is not consistent with a recent T2 measurement on YBa2Cu40s ||^. The Xs and 
Cu 1/TiT for Y Ba2CuiOs are similar to those of YBa2Cu30Q,Q but the T2 measurement implies 
that / dqlx'iq)]"^ ~ Xaf increases monotonically and smoothly by a factor of 2 between 300K and 
lOOK. If confirmed this would imply that as T decreases in YBa2Cu40s spin fluctuation weight 
is not only pushed away from low frequencies but also pulled down from high frequencies. Our 
simple model does not contain this physics. 

A realistic theory must incorporate itinerant carriers, and as in La2-xSr,j.Cu04^ a fermi-surface- 
instability description of the magnetic dynamics is required. One possibility, a model of two planes 
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of fermions, with direct hopping from plane to plane forbidden, was discussed in ref |T0[- Using a 
gauge theory formalism it was shown that an arbitrary weak interaction J2 eof Eq. 1 leads to a BCS 
instability to a "superconducting" state in which every Cooper pair has one member in each plane. 
We believe this superconducting state is a mathematical representation of the between-planes 
singlet; as pointed out in [jlOl it need not imply the presence of true superconducting order. We have 
not yet incorporated antiferromagnetism in the formalism of Ref. |T0|, but have shown that direct 
between-planes hopping of electrons acts as a pairbreaker of strength proportional to the charge 
carrier density [^. This provides a possible explanation of the difference between YBa2Cu30Q,Q 
and Y Ba2Cu^0^: in the latter material the much larger hole density permits interplane hopping 
which is strong enough to destroy the interplane pairing. 

In this paper we have proposed explanations for the spin gap behavior observed in underdoped 
cuprates. A crucial datum for our interpretation is the T ^ extrapolation of Xs{T)- We have 
argued that it is non-zero and indeed large in La2-xSrrcCu04 . If it is small, then a one-plane 
quantum disordered phase must be considered for La2-xSrxCu04, and the evidence for interplane 
pairing is in YBa2Cu30Q Q is weakened. Two important consequences are: (a) in La2-.xSrxCuOi 
samples with a T-dependent Xs the oxygen 1/TiT should decrease more rapidly than Xs with 
decreasing T and (b) in Y Ba2Cu^0Q,^ the Y relaxation rate should drop more rapidly than the O 
relaxation rate as T decreases. 
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Figure Captions 

Fig. 1.: Spin suspectiblities of liigli Tc materials obtained from data as described in the text, along 
with the theoretical susceptiblity of the 2d S=l/2 heisenberg model [|l8l. The powder average of 
the observed Xb for La2Cu04 at at T = 750K is within 0.2 states/ (ev-Cu) of the theoretical 
result. 



Fig. 2.: Copper NQR relaxation rates of high-Tc materials,from Refs. |T9|, |20|. Ref. pO[] used a 



normalization convention which differs from one used here by a factor of 3. We have reexpressed 
the data of ref p^j accordingly. 



Fig. 3.: Copper, oxygen and yttrium relaxation rates calculated for model of two coupled antifer- 
romagnetically correlated planes using Schwinger boson mean field analysis of eq 1 for J2 — J* = 0.3. 
The left ordinate shows the Cu and O relaxation rates 1/TiT (solid lines); the right ordinate shows 
the ratio of the O and Y 1/TiT to the calculated spin suceptibility Xs- 
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